Clustering and collision of inertial particles in random velocity fields 
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The influence of clustering on the collision rate of inertial particles in a smooth random velocity 
field, mimicking the smaller scales of a turbulent flow, is analyzed. For small values of the the 
ratio between the relaxation time of the particle velocity and the characteristic time of the field, 
the effect of clusters is to make more energetic collisions less likely. The result is independent of the 
flow dimensionality and is due only to the origin of collisions in the process of caustic formation. 
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The transport of finite size particles in turbulent flows 
is a common occurrence in several environments: rain- 
drops in clouds plankton in oceans 0, sprays in in- 
dustrial flows [3j, to make some examples. Due to inertia, 
these particles undergo clustering phenomena that have 
been observed in numerical simulations 0] , experiments 
and have been the subject of substantial theoretical 
study @, 0,1, 51 E|- 

Although contributing to particle segregation [l2| , spa- 
tial inhomogeneities in the turbulence do not appear to 
be an essential factor. What seems to be important is the 
ability of the particles to catch one another in their mo- 
tion, as they slip with respect to the fluid, a circumstance 
that is most evident in one dimension (ID) (HQ. In 
more than ID, the process is more complicated and an 
important role is played by the preferential concentration 
of heavy (light) particles in the strain (vortical) regions 
of the flow [6|. 

An important motivation for the interest in clustering 
is the contribution to particle collision and coalescence, 
and it has been suggested that this is an important in- 
gredient in the process of rain formation p], [ll| . What 
is observed is the simultaneous onset of concentration 
fluctuations and increased collision rates, when the re- 
laxation time of the particle velocity relative to the fluid 
becomes of the order of the turnover time of the fastest 
turbulent eddies [lj|. For sufficiently small (and suffi- 
ciently dense) spherical particles, this relaxation time is 
the Stokes time ts = 1/18 a 2 A/Vo, where a is the parti- 
cle diameter, A is the ratio of the particle to fluid density 
and vq is the kinematic viscosity of the fluid [l~6l |. 

When inertia is sufficiently high, the so called sling 
effect ensues [l7|: particle and fluid trajectories detach 
on the scale of eddies with turnover time rs, and their 
velocity will determine the particle collision velocity. 

The collision rate R co ii depends both on the particle 
concentration n(x, t) and on the relative velocity v in the 
particle pairs; for binary collisions: 

Rcoii ~ a 2 (n(a)n(O)) {v\r = a), (1) 

with r the particle separation. (Similar expressions will 
hold for the coalescence rate, with a more complicate 
function of u in the conditional average). Thus, both 
the sling effect and clustering, through the factors (v\v) 
and (n(a)n(O)), respectively, may be expected to enhance 



collisions. It was suggested in [Tt], EH, however, that the 
dominant contribution to collision may be the sling effect 
and there is even some indication (l9| that clustering may 
hinder rather than enhance collisions. 

The problem of how clustering and the dynamics of 
the two-particle velocity distribution in general influence 
each other, is far from trivial, with caustic formation po- 
tentially playing an important role [IT], [2(| • Purpose of 
this rapid communication is to understand whether it is 
possible to identify a clustering contribution to the ve- 
locity dynamics, and if this is associated with collision 
enhancement or hindering. 

Let us consider the dynamics of an inertial particle sus- 
pension in a smooth, incompressible random velocity field 
u(x, i), with correlation time te, variance (u^) = cr 2 , 
a = 1,2,3 and correlation length r v ~ <j u te- Particles in 
a turbulent flow, with rs shorter than the Kolmogorov 
time, here identified with te, will see turbulence precisely 
in this way, and this is appropriate for most aerosols of 
atmospheric interest [l|, [jj]]. Conversely, the opposite 
regime r§ > te could be interpreted as a model for the 
effect of eddies at scale r v on particles with ts corre- 
sponding to the turnover time of larger eddies. This limit 
is more relevant for industrial flows than for rain forma- 
tion, as the relative motion of larger droplets in clouds is 
dominated by the different gravitational settling veloci- 
ties of droplets of different size [l[ . 

The two regimes of large and small Stokes number 
S = ts/te are qualitatively differents, with the clus- 
tering maximum occurring somewhere at S < 1 [H, [2l[ . 
For small S, the particle phase is monodisperse in ve- 
locity over most of the fluid volume 0, [!]• The sling 
effect, due to the spatial correlation of u, occurs in co- 
herent way through the formation of caustics, i.e. re- 
gions of crossing of particle jets with different velocity 
[171] . The two-particle dynamics, relevant for the descrip- 
tion of clustering and binary collision, is described by 
the equation for the velocity difference T$i> + v = A r u, 
A r u(x, t) = u(x + r,<)-u(x, t), and, for {v 2 \r) ~ 

(|A r u| 2 ) -> for r -> (clearly, (z/ 2 |oo) = a 2 v = 2ct 2 ). 
Neglecting caustics would lead to O(S) collision rates, as 
predicted in the theory of Saffman and Turner [2^ . The 
multivalued velocity distribution in caustics appears to 
be crucial in producing high enough collision rates. 

In the opposite limit S> 1, the particles are scattered 
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by the velocity fluctuations they cross in their motion 
as if undergoing Brownian diffusion [23]. In the 5> 1 
limit, the velocity difference equation can be approxi- 
mated by a stochastic differential equation (SDE) in the 
form, choosing units such that t<j = a u = 1: 

v a + v a = b a p(r)£i3, r a = v a , (2) 

where 6 Q7 (r)6 7/3 (r) = (4/S)\g a p(0) - ^(r)], g a p(r) = 
a~ 2 (u a (r,t)up(0,t)), and £ Q is the white noise: 
(£a(*)£/8(0)) = S a p6(t). In the incompressible case: 
d a g a p(r) = 0. For S 3> 1, the majority of particles 
pairs at r = 0, will have been at r > r v for most of 
the previous time interval of length ~ t$ of which they 
have memory, so that it is possible to set, in the first 
approximation: b(r) ~ b(oo). This leads to a Brownian 
collision dynamics, with a velocity distribution of width 

Let us consider the effect of clusters on the collision 
velocity of the particle pairs, in the small Stokes numer 
regime. This regime could be analyzed within an SDE 
approach, imposing an artificially short correlation time 
to the random field te -C cr u /r v . This is a Kraichnan 
model regime [24j . in which the role of effective corre- 
lation time is played by the diffusion time for a pair of 
tracers, te — r 2 / (ct 2 te), to reach separation r v . We have 
therefore an effective Stokes number 

e = t s /t e = T S TE<rl/rl, (3) 

and it is possible to have a small-e large-S 1 regime, in 
which Eq. ([2]) continues to be valid. The small e 
regime has been the subject of extensive study (see e.g. 
[Til . [Til . |20T]). For small e, the particle phase is still 
monodisperse away from caustics, only with velocity not 
locally equal to that of the random field, as is instead (to 
first approximation) in the S <C 1, te ~ r v /a u regime. 

For small e, Eq. @ leads to r changing little in a 
time ts (the correlation time for u) and caustics arise 
as extreme events, in which a strong fluctuation in the 
random field causes particle pairs to jump ballistically 
to zero separation in a time ~ rg fl4l |20T ] . Clusters af- 
fect the process privileging particle pairs that are initially 
closer, i.e. less energetic fluctuations in A r u, leading to 
smaller collision velocities [25| . We give a quantitative 
description of this effect for D = 1 and |r| <C r v . 

For e — > 0, the shapes of the pair trajectories termi- 
nating with a given collision velocity, will concentrate 
around the one that maximizes probability (no condition 
is imposed on the caustics in which the trajectories de- 
velop). To determine the most likely trajectory ending 
with collision velocity v(t) — v at t = 0, one can proceed 
iteratively, from some initial guess for the separation his- 
tory, say ro(t) = vt. (We assume immaterial particles, so 
that they can overlap without interaction) . For a smooth 
field, we can take for |r| <C r v g(r) ~ 1 — a(r/r v ) 2 with 
a = 0(1), which gives b(r) ~ e 1 / 2 r. For r -c r v , at 
the fc-th step in the iteration procedure, the first of Eq. 



([2]) will then read, apart of an O(l) factor in front of the 
right hand side (RHS): v k + v k = e^ 2 ?"*^, leading to the 
solution 

Mt) = t 1/2 I dre T -V fe (r)^(r), (4) 

J — oc 

where = Vk-i, ^fe(O) = 0, k > 0. The minimum 
problem will be, therefore, 

5pog^K] + A k i/ fc (0)] = Q, (5) 

where S indicates variation in £, V[£\ is the functional 
PDF of the history £(t) with t £ [— oo, 0], Xk is the La- 
grange multiplier to enforce ffe(0) = v and the entering 
Vu is assigned from the previous iteration. From V[£] — 
cxp(— i J £ 2 (t)dt) and Eq. ([U, the minimum problem 
Eq. © leads to the history £ k (t) = \ k e 1/2 e* r k (t). 

Substituting into Eq. ((4]) and imposing 2/(0) = v 
gives the velocity profile v\~{t) in function of r k , which, 
substituting into rfc+i = v k leads to the iterative relation 
for r k (t): 

r k+1 {t) =9{l-[J oo dre 2 W k (r)}^ 

The iterative scheme can be implemented numerically 
and converges rapidly. In particular, the starting point 
of a jump terminating in a collision at velocity v is r p ~ 
linifc^oo r k (— oo) and one finds from Eq. (JSJ) r p ~ —Ivrs- 
The rather large factor 7 implies that the jump does not 
start from a spatially localized "kick", rather, the sling 
acts over an 0(Dts) distance equivalent to that of the 
final free flight. 

From the same condition is(0) — one finds X k — 
[e J_ dr e 2T r k (T)]~ 1 D, which, substituted into £fc(i) = 

Afee 1 / 2 e* r k (t), together with r k ~ D, gives for the noise 
history: ~ e" 1 / 2 and therefore for the PDF of ob- 
serving a collision velocity v in a given pair will be, for 
a < r p : Pj ump {D\a) ~ p 3ump (v\0) ~ ~ exp(-c/e), 

with c a constant that can be shown to be 1 /6 in ID [lj] . 
The PDF of the collision velocity v in the given pair, 
generated by a jump originating at separation r will be 
instead: Pj U mp (v\0)d(r — vT) where T ~ 7ts- 

Multiplying by the PDF p(r) of finding a pair at sep- 
aration r, and integrating over r, gives the PDF of a 
collision at velocity v 

Pcoii(v) = p(r^)pj um p, (7) 

with the cluster contribution contained in the PDF for 
the particle pair separation p(rp). A collision velocity v 
implies a permanence time oc \v\ _1 in an interval dr, and 
this allows to write the collision velocity PDF (that is 
the collision rate at that velocity) in the form 

Pcou{v) - p(D, a)\9\ = p(a)p{v\a)\v\. (8) 

This in turn can be substituted into Eq. (TT]) exploiting 
the relation, valid in D generic: (n(r)n(O)) = fln 2 p(r), 
with £1 the domain volume and n the mean concentration. 
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FIG. 1: Collision velocity PDF for the ID problem u + u = 
e 1//2 (?- + a)£, r = u, compared to p = v~ 2 (upper line). Values 
of the parameters: e = 0.04 and a = 0.02 (a finite molecular 
diffusivity is necessary to regularize the dynamics at r — * 0). 



In ID, p(r) oc r~ 2 pjj], so that p co uiv) oc v~ 2 , as 
confirmed in Fig[TJ In the absence of clustering, in com- 
parison, the distribution would have been uniform and 
~ exp(— c/e) up to \v\ ~ r v /r s . 

This picture extends to D > 1, as the extremal trajec- 
tories are still straight lines, and the relevant parameter 
remains the separation r p at the start of the jump. In 
this case the finite particle size must be taken into ac- 
count and a jump ending in a collision will develop along 
a straight line that does not necessarily pass through the 
center of the other particle. The probability of a jump 
originating at r leading to collision will be proportional 
therefore to the angle with which a particle is seen at 
distance r, i.e. (a/r) ^ 1 . Nevertheless, if a -C rp, the 
collision velocity PDF, i.e. the PDF for the velocity in 
the direction of the jump at the other particle position, 
could be approximated as: Pj um p(v\a) ~ pjump{ v \^)- 

Now, the number of particles in a shell at distance 
r will be oc p{r)r D ~ l dr, where the PDF p(r), pro- 
vided lim r ^o r D p(r) < oo, is related to the correla- 
tion dimension Di of the distribution by the equation 
p(r)r D oc r~ D (N 2 ) oc r° 2 , with N r is the number of par- 
ticles in a volume of linear size r. Taking the product of 
the different contributions, the probability of collisions 
at velocity between v and v + dv will be therefore 

p C oii(v)dv oc p ]ump {D\Q)p(r 9 )a D ~ 1 dr i)l (9) 

where r p = vT for some T = 0(rg). If, as in the ID 
case, Pjump(v\0) is independent of is, the following result 
will hold: 

Pcoii{v) « D° 2 - D . (10) 

Simulating the trajectories of an ensemble of inertial par- 
ticles in a 2D Kraichnan random field leads to the result 
in Fig. [21 which confirms the prediction of Eq. (fit)]) . 
The peak to the left is produced by the finite size of the 
particles, in the present case a ~ 1/2056 of the domain 




FIG. 2: Collision velocity PDF in 2D from a sample of 10 5 
particles with e = 0.04, advected by a smooth Kraichnan ran- 
dom field in a square domain 3r v x 3r v . The almost straight 
line is p(r) = r~ 2 (N 2 ) (appropriately rescaled), correspond- 
ing to a correlation dimension D2 — 1.45. 




FIG. 3: A snapshot of the particle distribution from the 
numerical simulation of Fig. [2] superinposed to the pattern 
of relative velocities (the small segments) , which identify the 
caustics. The length and direction of the segments correspond 
to the largest eigenvalue and associated eigenvector of the 
correlation matrix (i/iVj\x.,t), calculated as average on the 
particles in a box of size a at x. 



size. Its width is (z^la} 1 / 2 ~ aa u /r v (the typical rela- 
tive velocity at separation o) and its height is a factor 
~ exp(c/e) above the scaling range to the right, which is 
associated with the jumps. 

As illustrated in Fig. [3l the collisions actually take 
place between clusters. The scaling D2 < D causes 
closer clusters and therefore less energetic collisions to 
be more likely. Notice that, although clustering hin- 
ders high velocity collisions, higher velocities are actu- 
ally more probable inside clusters. The reason is purely 
statistical: higher particle concentrations are produced 
where clusters collide, i.e. by definition the place where 
v is larger. 

For large Stokes numbers, collisions cease to occur as 
extreme events; hence, the correspondence between colli- 
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sion velocities and jump lengths disappears, and the pre- 
dictions of Eqs. (|7I10[) cease to be valid. In this limit, the 
particle velocity distribution ceases to be monodisperse 
pointwise (near caustics, it would be a superposition of 
discrete jets) and particle velocities at close separations 
are in the first approximation independent. It was sug- 
gested in flU] that concentration fluctuations are pro- 
duced by slowly approaching particle pairs, which spend 
a significant time at separations r < r v . To verify this, 
however, it is necessary to separate out a cluster contri- 
bution in the velocity PDF p(v\r) 1 analogous to the one 
identified in Eqs. (fST TTUj) . or equivalently in the relation 
p{D,a) oc 1 1> | 1 p um P {v\a)p(r 9 ). 

Notice that the simultaneous increase of p(r) and de- 
crease in (v 2 \r) as r — > is not sufficient to conclude that 
concentration fluctuations are associated with smaller 



relative velocities; in fact, lim r _ > o( I/2 | r ) = also for pas- 
sive scalars in incompressible flows, in which case, con- 
centration fluctuations are absent. An alternative ap- 
proach is therefore required. 

As carried on in [261 ] . a clustering part in p(y,r) — 
(8(u(t) — p)8(r(t) — r)} could be identified in the higher 
order contributions from the expansion of r (t) a round the 
Brownian motion limit. The analysis in [26J] indicates 
a strong dependence on compressibility of the random 
field, and that the conclusion of the present paper, that 
clustering decreases collision velocities, remains valid at 
large S only for compressible flows. As in the small e 
regime, however, this clustering contribution cannot be 
identified with the actual velocity distribution inside the 
clusters, and is more in the form of a non-local cluster 
contribution to the velocity dynamics. 
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